Single-photon, cavity-mediated gates: detuning, losses, and non-adiabatic effects 
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We study several extensions of the single-photon, cavity-mediated quantum logical gates recently 
proposed by Koshino, Ishizaka and Nakamura: to a double-sided cavity configuration, to the case 
where the two atomic ground states are nondegenerate, and to the non-adiabatic regime. Our 
analysis can be used to estimate the effects of various imperfections, and to prepare the way for a 
proof-of-principle demonstration with present technology. We are able to present a full solution for 
the outgoing pulses, in the frequency domain, which is valid for arbitrary incident pulses and cavity 
parameters. From this we find, among other results, that the leading correction to the adiabatic 
approximation can be made to vanish for a suitable choice of detunings, provided the cavity is "good 
enough" (high enough ratio of coupling to loss). This could significantly relax the need for long 
single-photon pulses in this scheme. 



I. INTRODUCTION 

Photons are widely regarded as natural carriers for 
quantum information, but their usefulness is limited by 
the lack of strong photon-photon interactions that could 
be used to carry out conditional quantum logic on single- 
photon states. Although Kerr nonlinearities have often 
been proposed to address this difficulty, there are good 
reasons to believe that they cannot really be useful for 
quantum logic at the single photon level [l|, , at least 
in an unconstrained geometry Q. 

Atomic ensembles and optical cavities can be used to 
constrain the interaction geometry and increase the ef- 
fective photon-matter coupling, thus providing a possi- 
ble way to mediate single-photon interactions. A single- 
photon gate mediated by a single, three-level atom in a 
cavity was proposed by Duan and Kimble in 2004 4] and 
generated a great deal of interest, but it has not been ex- 
perimentally realized yet because of its very challenging 
nature. Part of this challenge is the need to address the 
atom inside the high-finesse optical microcavity with an 
external laser field (presumably through the sides of the 
cavity), in order to perform the gate. 

In general, all proposed variations of the Duan-Kimble 
scheme rely on the different phase shifts acquired by a 
photon upon reflection from the cavity, depending on 
whether it is resonant or not with the atomic transi- 
tion. Although, as indicated above, a full gate of the 
Duan-Kimble type has not been demonstrated yet, the 
required conditional phase shifts have by now been ob- 
served in a variety of systems Q , and many potential 
uses for them in quantum information processing have 
been suggested (see, for a very small sample, |7H10|, and 
references therein). 

In contrast to the above proposals, it was shown in 
2009 by Lin et al. [llj that, if a A- type atom with two 
degenerate transitions (corresponding to different polar- 
izations) in a doubly resonant cavity was used, the inci- 
dent photon could directly change the atomic state, so 



that a SWAP gate between the atom and the photon po- 
larization states was possible. More recently, Koshino, 
Ishizaka and Nakamura (KIN) showed [12| that a totally 
passive, cavity-mediated gate between photons was pos- 
sible in the same system for an appropriate choice of de- 
tunings. In their scheme, a VSWAP gate between two 
single-photon pulses (a two-qubit gate which, together 
with single qubit gates, suffices to enable universal quan- 
tum computation [l3j ) could be carried out merely by 
reflecting the pulses off of the cavity in an appropriate 
sequence: manipulation of the atom by external fields, 
whether lasers or static, is not required at all, in princi- 
ple. Moreover (as already proposed in [11]), the KIN 
system could also be used — again in a totally passive 
way — as a quantum memory, by swapping the photonic 
and atomic states. KIN modules could therefore pro- 
vide a useful alternative to atomic ensembles, both for 
storage and manipulation of quantum information. In a 
previous paper [l4j we have showed how the KIN scheme 
(originally proposed for the "bad" or "fast" cavity case) 
could be extended to the "good" cavity regime. Other 
properties of the KIN system, in particular its ability to 
function as an ideal photon turnstile, have also been re- 
cently reported [Hj]. 

In this paper we wish to extend our investigation of 
the KIN system by relaxing some of the assumptions 
and restrictions on which previous studies (including our 
own) have been based, with an aim towards ascertaining 
the feasibility of a proof-of-principle demonstration with 
present technology. The outline of the paper is as follows. 

In Section II, for completeness, we present a very brief 
overview of the original KIN proposal. In Section III, 
we consider what happens if one replaces the original 
single-sided cavity by a double-sided cavity, an impor- 
tant question given that many of the best current opti- 
cal microcavities are either double-sided or ring cavities 
(such as microdisks, or microtoroids), and we would like 
to explore the possibility of using an existing cavity for 
at least a proof-of-principle experiment (note that some 



previous proposals, such as 10], naturally work in the 
double-sided geometry, so it is perhaps not immediately 
obvious why the KIN scheme would not). Formally, also, 
a "leaky" second mirror can be used to model the effect 
of cavity losses other than those arising from the trans- 
mission through the first mirror. 

In Section IV, we explore the case in which the atomic 
transitions are not degenerate. This too is important for 
several reasons: for example, in some atomic systems the 
only way to select a closed, lambda-like set of three hy- 
perfine states may be by using external magnetic helds to 
make all unwanted transitions non-resonant with the fre- 
quencies of the photons used. In such a scheme, the two 
polarizations making up the photonic qubit would nat- 
urally have to have different frequencies as well (as we 
show in detail below). Again, alternatively, our results 
may be taken to show the damaging effect of atomic en- 
ergy shifts if these are not compensated for. 

Lastly, in Section V, we consider the corrections to the 
adiabatic approximation under which all the KIN results 
so far have been obtained. Interestingly, we show that the 
leading correction to the adiabatic approximation can be 
made to vanish for a suitable choice of cavity parameters, 
which could significantly relax the need for long pulses. 



II. THE ORIGINAL KIN SCHEME 

In the original KIN proposal, one considers a single 
three-level atom, in the A configuration, inside a single- 
sided cavity (Fig. 1). The two ground states are taken to 
be degenerate, so the two transitions shown correspond 
to different polarizations. A single-photon pulse enters 
the cavity, interacts with the atom, and leaves through 
the same input mirror. In the process, it may change 
the state of the atom or become entangled with it. If 
the cavity decay rate is n, the atomic decay rate is 7, 
and the single-photon Rabi frequency (coupling constant 
to the cavity mode, assumed to be the same for both 
transitions) is g, the "strong coupling regime" (assumed 
throughout) is that 7 <C n,g. Within this regime one 
may still distinguish a "bad-cavity" limit, k>j, and a 
"good cavity limit," n <C g. In the bad-cavity limit, the 
cavity parameters enter the equations only through the 
combination 2g 2 /n, and the strong coupling condition 
becomes 7 <C 2g 2 / n. 

In the bad-cavity limit, Koshino, Ishizaka and Naka- 
mura (KIN) showed [12| that, for a sufficiently long pulse 
(adiabatic regime), the input and output states of the 
total atom+photon system were related by the transfor- 
mation 




\H,0)->\H,Q) (la) 

\H,1) -> -e^ [iaisi(f>\H,l) + cos(f>\V,0)] (lb) 

|V, 0) -> -e^[cos0|H,l) +isin0|KO)] (lc) 

\V, 1) -> \V, 1) (Id) 




FIG. 1: The KIN scheme. Single-photon pulses are reflected 
off of a one-sided cavity containing a three-level atom, in the 
process becoming entangled with the atom. By swapping the 
state of a second pulse with that of the atom, an entangling 
operation between single-photon pulses results. 



where is a detuning-dependent parameter. When 
= one has essentially (up to a sign) a SWAP gate be- 
tween the photon and the atom, whereas when cf> = 7r/4 
one has a VSWAP gate [lj]. To have <j) = it is suf- 
ficient that all the detunings vanish. To have <fi — 7r/4, 
Koshino, Ishizaka and Nakamura derived the condition 
5 a = 2g 2 / k for the detuning between the atom and the 
cavity, assuming implicitly that the cavity and the field 
were resonant. For the more general case where the field 
is detuned from the cavity resonance by an amount A, 
we showed in [l4| that the condition to have (f> — it/ 4 is 



(A + <5 a )(« 2 
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This also reduces to the KIN condition in the very bad 
cavity limit, k» A. The above results will be seen to 
follow as special cases of the more general analysis to be 
carried out later in this paper. 

In their original work, Koshino, Ishizaka and Naka- 
mura studied the effect of some deviations from their 
assumed ideal setup, such as different coupling constants 
for the horizontal and vertical polarizations, and finite 
pulse length. In as stated above, we extended their 
results to the good cavity limit, considered (in a limited 
way) the effect of spontaneous emission, and established 
(numerically) the 1/T 2 dependence of the gate fidelity on 
pulse duration T. 

In the sections that follow, we shall study, mostly an- 
alytically, the effects of other departures from the ideal 
KIN assumptions: two-sided cavity (or additional cavity 
losses) ; nondegenerate ground state; and deviations from 
the adiabatic evolution (finite-pulse effects). Our work 
will be quite self-contained, since we will even present an 
alternative derivation, and extension, of our "modes of 
the universe" formalism in the following section. 



III. KIN GATES WITH A DOUBLE-SIDED 
CAVITY 

The "modes of the universe" approach was originally 
introduced by Lang, Scully and Lamb [l6| as a rigorous 
way to quantize the quasi-modes of a lossycavity, such 
as a laser cavity. The formal treatment in [171 , on which 
our previous paper was based, involves a number of tech- 
nical subtleties and is not physically very intuitive, so in 
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FIG. 2: The scattering modes of a double-sided cavity 



for all the relevant modes, 



2ikl 



where Q c is the cavity's resonant frequency, in the de- 
nominator of whereas in the numerator we can sim- 
ply replace e lkl by —1. The reason for this is that, for 
a microcavity, c/l is of the order of an optical frequency, 
and the typical detunings we shall consider will be smaller 
by many orders of magnitude. 

We can also assume that the reflectivity of the mir- 
rors is very close to 1 and that they are lossless, so 
r = Vl - t 2 ~ 1 - t 2 /2. Substituting all this in Eq.® 
we get for the intracavity field 



what follows we present an alternative derivation which 
is also an extension, to the case of a two-sided cavity 
with mirrors of (amplitude) transmission coefficients t\ 
and t2- 

The idea is to consider the coupling of an atom inside 
the cavity to the "scattering modes" depicted in Figure 1. 
For a given frequency fife = ck there are two such modes, 
one incident from the left and one from the right. For 
unit incident amplitude, the resulting field amplitudes 
illustrated in Fig. 2 are, for the mode incident from the 
left, 
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r\r2e 2lkl — 1 



(3) 

assuming the cavity has length / and the mirrors have 
real reflection coefficients, whereas for the mode incident 
from the right the coefficients A 1 , B' , C and D 1 are given 
by identical expressions with the indices 1 and 2 switched 
around. 

Let a>k (ct_fc) be the annihilation operator for the mode 
incident from the left (right). At an arbitrary point inside 
the cavity, an atom would see a field given by 



\ - / hrt k 



1/2 
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_L + f„ J-kl ikx —ikx\ 

+ £2 me e — e I a_fc 



(4) 



where AL is the quantization volume for the scattering 
modes. The resonance condition kl = (2n + 1)tt mini- 
mizes the denominator in Q while keeping the numera- 
tor close to a maximum (for an atom near the center of 
the cavity, x ~ 0). We assume then that we can write, 
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k - i(p,k - n c ) 



(6) 



to lowest order in small quantities. In the last equality 
we have replaced Q k by the pulse's central frequency Qo 
in the "electric field per photon" factor. We have also 
introduced the total cavity loss rate 



tf + q c 
4 7 



and the parameters 



+ t 



1,2) 



(7) 



(8) 



satisfying + r| = 1. 

Suppose one has a two-level atom at x = 0, and let 
the coupling energy be of the form Ed (where d is the 
atomic dipole moment); then, in the rotating- wave ap- 
proximation, the interaction Hamiltonian can be writ- 



ten as Tig{a cajV (j* 



where the coupling constant 



g = (d 2 Qo/heoAl) 1 / 2 , is an atomic raising operator, 
and the intracavity effective quasi-modc field operator is 



E 



v / 2kc7X 



k - i(Qk - fi c ) 



(nafc + r 2 a_ fc ) (9) 



One can verify that [a c 



J = 1, using the canon- 



ical commutation relations for the a k ,a^k modes, and 
converting the sum over k to an integral via J^ fe — > 
2ir/L J dk; here 2irc/L is the frequency spacing appro- 
priate for running waves in a cavity of length L (with 
periodic boundary conditions), as opposed to irc/L for 
the standing- wave quantization assumed in [lij . 

The crucial point of Eq. © is that, in general, the field 
that couples to the atom is a superposition of modes in- 
cident on the cavity from both sides. Defining "coupled 
mode operators" a ck — T\a k + T2<2_fc and "uncoupled 
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mode operators" a u k = r^a^ — ria_fe, it is clear that 
if the cavity is driven only from, say, the left side by a 
one-photon pulse, the incident state is a superposition of 
the form ti|1) c + t 2 |1)„, where the second term does not 
participate in the interaction at all. If our goal is, for in- 
stance, to change the state of the atom (as with a SWAP 
gate) this immediately leads to a failure probability of at 
least r|. 

Essentially, "opening" the other side of the cavity turns 
a scheme for quantum logic that was originally deter- 
ministic into a probabilistic one. The damage can be 
quantified in the following way: in Eq. (TTJ, take all the 
photon states on the left-hand side to be left-incident 
single-photon states, and write them as a superposition 
of coupled and uncoupled modes, as above; apply the 
transformation (TTJ) to the coupled part only; then rewrite 
the final result in terms of left- and right-incident modes 
again. The "error" terms are all those containing a factor 
of r 2 . 

These results highlight an important difference be- 
tween the KIN scheme and the Duan-Kimble-type gates. 
The latter always work in the limit of negligible atomic 
excitation, in which the atom-cavity system acts as a 
whole, imparting a phase change to the photon without 
changing its own state. In the KIN scheme, instead, it is 
essential that the photon couple strongly enough to the 
atom to change its state, and in general it is impossi- 
ble for a single photon coming from a definite direction 
to deterministically change the state of an atom that is 
coupled just as strongly to a vacuum mode traveling in 
the opposite direction. The same problem would be en- 
countered, therefore, in ring cavities such as microdisks 
or microtoroids. 

In principle, however, as long as the cavity is lossless 
a way around the difficulty exists. One must, in the first 
place, drive the cavity from both ends in such a way as 
to excite only the coupled mode, and then collect all the 
light leaving the cavity and recombine it appropriately. 

A possible arrangement is shown in Figure 3. For op- 
timal coupling, the beam splitter should transmit mode 
a' with amplitude n, and reflect it with amplitude r 2 . 
Then, as shown in the figure, one has = T\a! — T 2 a" , 
and a_fe = T±a" + r 2 a' (where a" is a mode in the 
vacuum state), and the "coupling mode" combination 
&ck = Tidk + T2a.-k = a'- Hence a photon entering the 
beam splitter at port a' goes whole into the coupling 
mode. 

The atom also emits the photon in the coupling mode, 
so its fate upon reaching the beamsplitter can be de- 
termined by considering the linear superposition of the 
four waves A, D, A' and D' in Figure 2, with the pre- 
cise phase relationship corresponding to that mode. Note 
that in our solution we have assumed, as shown in Fig. 1, 
that the incident field at the right mirror and the incident 
field at the left mirror have the same phase, so in Figure 
2 with Ti and r 2 real and positive the optical path lengths 
from the beamsplitter to either mirror must be the same. 
Then the left-traveling field arriving at the beam splitter 



a ' > r x 
— *\ 



V '-\ a 'j ( — \ 

> — f v/\h| 

BS1 N l°> I 1 ) 7 
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FIG. 3: (Color online) Optimal coupling to the double-sided 
cavity (see text for details). The lighter arrows represent the 
outgoing field. The incoming field a" is in the vacuum state. 



by the upper path is proportional to T1A + T2D' , and the 
one arriving from below is proportional to t\D + t%A! , 
times a common phase factor. Near resonance we have 
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(10) 



Using these expressions, we find that the left-traveling 
and upward-traveling fields at the beamsplitter go as t\ 
and T2, respectively, times a common phase factor. Com- 
bining them, we see that no field leaves the beamsplitter 
in the upwards direction (since the reflection coefficient 
on that side of the mirror must be equal to — t%). Hence, 
in this arrangement, a photon emitted in the coupling 
mode must leave the beam-splitter along the same path 
(but in the opposite direction) as the input field. 

We conclude that the arrangement in Fig. 2 allows 
one to couple optimally to the cavity-atom system, for 
both input and output. A similar arrangement would be 
possible for a ring cavity, for which the coupling mode 
would be a standing wave; as in Fig. 2, in that case also 
one should use a beamsplitter to feed the input field into 
both counterpropagating modes simultaneously, and to 
collect the output field. 



IV. NON-DEGENERATE GROUND STATE 

Returning to Eq. (j9|), suppose that either the cavity 
is one-sided, so r 2 = 0, or we arrange to excite, and 
collect, only the coupled mode. In either case we can 
then just use the symbol ak for that mode's component at 
frequency f2fc = ck. Introducing the two polarizations H 
and V, and an atom with two corresponding transitions, 
we find, in the interaction picture, the Hamiltonian for 
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the system: 



y/2cn/L 



k - i(Q k - n c ) 



(a hfc |2)(l|e 



—i(Q k —u! h)t 



+ a„ fc |2)(0|e^ (o ^"°" )t ) + H.c. (11) 

where the atomic states are labeled as in figures 1 and 
3 and u)Qh an d u>ov are the corresponding resonant fre- 
quencies. This Hamiltonian differs from the one we 
used in [l4j in only a few particulars: the quantization 
length/frequency spacing (already mentioned above), a 
phase factor which is discussed in more detail in the Ap- 
pendix, and the possibility of nondegenerate atomic tran- 
sitions, with which this section is concerned. 

The original KIN proposal assumes that the two 
ground-state sub levels that make up the atomic qubit 
are degenerate. This may not always be the case; for 
instance, for multilevel atoms one might want to delib- 
erately break up the degeneracy of the various magnetic 
sublevels in order to confine the evolution to an appropri- 
ate manifold. Other systems that have been proposed for 
quantum logic in a solid-state setting explicitly start out 
with two non-degenerate lower levels: for instance, dou- 
bly excited quantum dots, or nitrogen vacancy centers in 
diamond. 

Assuming the state vector is of the form 



We shall begin by assuming that the pulse is horizontally 
polarized and the atom in the state 1 1) . We can integrate 
formally Eqs. (|13a|) and (|13bj) . substitute in (|13c|) . and 
carry out the sum over k in the continuum limit, as we 
did in [14j . The result is the integral equation 



C e = -g 2 [ (e-( K +^)( t -*')-|- e -( K + w ")(*- t '))c e (t')^' 



ig /d«V^(7 h ( Wj 0) e - i ("+^ 



(14) 

where the atom-cavity detunings Sh — s7 c — i^oh an d 
5 V = Q c — ujq v have been introduced, as well as the 
continuous version of the mode coefficients: Ch((jJ,t) — 
Chk(t)\/ L/2-kc (this is appropriate for a mode spacing 
2irc/L and ensures / \C h {uj)\ 2 duj = £ fc |C, lfc | 2 = 1). The 
frequency U) is taken to be zero at the cavity resonance. 
The advantage of the continuous formulation is that it 
leads to a straightforward solution of (IT4l) (something 
which we failed to appreciate in [HI]); in the frequency 
domain. Introducing the Fourier transform of C e (t) by 



C e (t) = j- I Cei^e-^du; 



(15) 



!*(*)> = ^Z G hk{t)\l)a J< hk \v&c) +^2c vk (t)\0)al k \va,c) it i s easy to sec that Eq. (O implies 



+ C e (t)|2)|vac) 



(12) 



C e (w) 



-2ig 



(where |vac) is the field vacuum state) we obtain the 
equations of motion 



C h {u-8 h ,0) 



g 2 (2 + ie/(n — i(u + S v ))) — iuj(k — i(u — Sh)) 



C h k = -ig 
Cvk = -ig 



k + i(Qk - sl c ) 

V 2cK / L e i(n k - Uov )t 
k + i(Ofc — sl c ) 



C e (t) 



(13a) 



C e (t) (13b) 



^ k — l(Mk — s2 c ) v 



(13c) 



(16) 



where e = Sh — S v . Also, the time integrals of the right- 
hand sides of Eqs. (|13a[) and (113b[) . from well before the 
interaction starts (which may alternatively be considered 
the time 0, or — oo) to well after it is over (formally 
plus infinity) are then just proportional to C e (uj + Sh) 
and C e (u! + <$„), respectively. Putting this together with 
Eq. (ITcH) immediately yields the frequency spectra of the 
outgoing horizontally and vertically polarized pulses: 



2g 2 K C h (uj,0) 

k + iuj g 2 (2 + ie/(K — i(u) + e))) — i(u) + Sh){n — ioS) 

2g 2 K C h {oJ-e,0) 

K + iu; g 2 (2 + ie/(n — iuj)) — i(u> + S v )(n — i(ui — e)) 

I 

These expressions show that the outgoing vertical pulse is centered at a frequency different from the one of the 



C h (w,oo) =C h (u,0) — , t — ;— ^ - — r - FT7 — - (17a) 

C v (w, oo) =- - ■ — 2 ■ tttt^ — ^ 77 , , ; — — (17b) 
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incoming horizontal pulse, as one would expect from con- 
servation of energy (see Fig. 4). 

Leaving the question of this frequency shift aside for 
a moment, we may consider what is required in order to 
get the ideal transformation indicated in Eqs. ((T|). Ide- 
ally, the final pulse should be a perfect replica of the ini- 
tial one, which would only be possible if the coefficients 
multiplying Ch(w,0) in Eqs. (fTT)) did not depend on u>. 
This is clearly not the case, but for a very long pulse, 
centered around some frequency the spectrum may 
be so sharp that it is a good approximation to evaluate 
those factors only at that frequency. This is essentially 
the adiabatic approximation, and departures from it will 
be considered in the next section. 

For the remainder of this section, therefore, we will 
treat the prefactors in Eqs. (fTT)) as constants evaluated 
at the frequency where C/j(a;,0) is maximum. Note that 
u) in these expressions is defined relative to the cavity 
frequency, so if the incoming pulse peaks at fife, the cor- 
responding value of uj in Eq. (|17al) is VLh — fi c = A^; also, 
in Eq. (|17b|) . Ch{u — e, 0) peaks at u = A/, + £, so it 
is at this frequency that the coefficient should be evalu- 
ated. Under these conditions, therefore, Eqs. (fTT]) yield 
the approximate transformation 



\H,1) -> a\H,l) + 0\V,O) 



(18) 



with 



2g k 

a = l — x 

k + iAh 



g 2 (2 + ie/(K - i(A h + e))) - i(A h + S h )( K - iA h ) 
K + iAh 



K + i(A h +e) 



(I -a) 



(19) 



It is now an easy matter to see that, in the degenerate 
case (e = 0), the condition @ does indeed produce a 
VSWAP gate, that is, Eq. (fib]), with = tt/4. On the 
other hand, in the nondegenerate case one can still make 
the coefficient a equal to the desired value, (1 — i)/2, 
provided the detunings are chosen to satisfy 



Ah = -2 



2g 2 



(20a) 
(20b) 



K {l + e 2 /4n 2 ) 

(note this solution is unique). In that case one also has 
ie/2 l + i 



l+i 



K + ie/2 2 



(21) 



which is also equal to the desired value, except for a phase 
factor e , with 8 = tan _1 (e//c). As we shall see below, 
this does not necessarily spoil the %/SWAP gate opera- 
tion. 

We still need to check the third of Eqs. (fTJ (the first and 
fourth are, of course, trivial). For this the initial state 
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FIG. 4: Detunings necessary for a %/SWAP gate. The tran- 
sition |0) — > 1 2} couples to vertically-polarized radiation with 
frequency Q v , and the transition |1) — > |2) to horizontally- 
polarized radiation with frequency Qh- The resonant fre- 
quency of the cavity is Q c . 



should be chosen to be | V, 0) , and, for consistency, the 
vertically-polarized pulse should be chosen to be centered 
at f2„ = flh + e, as indicated above. In that case, all the 
math leading to Eqs. (fT7|) still holds, only with v and h 
interchanged, and the sign of e reversed. Now the choice 
Ah = —e/2 implies A„ = O/j + e — f2 c = Ah + e = e/2, 
and we obtain 



IV, 0) 



1 



IV o) 



l + i 



HA) 



(22) 



which is again the desired result except for a phase fac- 
tor. The form of this factor shows that, indeed, the op- 
eration is unitarily equivalent to a VSWAP, since it is 
the opposite of the one appearing in front of | V, 0) in the 
transformation ([15]) (according to Eq. (fUJ). In practice, 
one could eliminate these factors by multiplying the H 
component of the incoming pulse by e lB before the inter- 
action, and by e~ %e afterwards. 

In short: the VSWAP gate can be performed (in the 
adiabatic approximation) exactly provided the vertically 
and horizontally polarized pulses have, from the start, 
different frequencies. This frequency difference is then 
preserved throughout the interaction. The precise de- 
tuning condition needed is illustrated in Fig. 4. 

While this sounds good, the problem with having dif- 
ferent frequencies for the two polarizations is that single- 
qubit gates — which in the degenerate case were simple 
polarization rotations — now become difficult to imple- 
ment. A 45 degree rotation applied to both fields, for 
instance, creates a superposition of fields with different 
frequencies at each polarization component. So for a 
single-qubit gate one would have to first separate the po- 
larizations, frequency-shift them so that they both have 
the same frequency, combine them and rotate them, sep- 
arate them again and frequency-shift them back to what 
they were, and finally combine them again. 

Alternatively, as long as one is only doing single-qubit 
operations, one could just keep the two polarizations de- 
generate, and only up- or down-shift them before they go 
into the cavity, and after they come out. We note that 
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frequency shifting of single photons that preserves po- 
larization entanglement has been recently demonstrated 
experimentally [181 ] . so this is definitely feasible, in prin- 
ciple. 

For small e it may be tempting to ignore the frequency 
shifts altogether and work with (nominally) degenerate 
polarizations anyway. In this case, there is a dual price to 
pay: a mis-overlap of the pulses in frequency space after 
each interaction with the cavity, which can be estimated 
as ~ (eT) 2 for a pulse of duration T, and an error in the 
gate itself. Both of these can be calculated together from 
the exact Eqs. (IT7|) . although the gate error may also be 
estimated from the approximate Eqs. ([T8j) and fp~9|) . For 
instance, suppose one simply tunes both f2^ and fl v to 
the cavity resonance Sl c , and sets the latter so that Sh and 
5 V have their correct values, as in Fig. 4 and Eq. (|20b)l . 
Then the error in the coefficients a and /3 goes as 



1 

An 



(23) 



to first order in e. Interestingly, this first-order error 
can be made to vanish by using a cavity that satisfies 
g = k/-\/2. Nonetheless, even under these conditions, 
one still finds a second-order gate error of the order of 
(e/«) 2 . Recalling that the adiabatic condition requires 
kT 3> 1, we see that the gate error will typically be 
smaller than the frequency mis-overlap error; that is, the 
basic criterion for the frequency split e to be negligible is 



e 2 T 2 



< 1 



(24) 



which may be difficult to satisfy in the adiabatic regime. 
The requirement makes sense, physically, since in this 
regime the pulse duration provides the smallest band- 
width in the system, hence the pulse itself (not the cav- 
ity) is the sharpest frequency discriminator. 



CORRECTIONS TO THE ADIABATIC 
APPROXIMATION 



From all of the above we see that it is important to try 
to determine how long the pulse has to be, in practice, 
before the results obtained under the adiabatic approx- 
imation cease to be valid. Put otherwise, we want to 
calculate the leading corrections to the "very long pulse 
approximations" that were made in [14j . 

This is actually not difficult, since we have the exact 
solution (|17l) in the frequency domain, valid for any initial 
pulse. For simplicity, in the remainder of this section we 
will be concerned only with the degenerate case, e = 0, 
and we will simply write 5 a for the atom-cavity detuning 



fl c - w - Then Eqs. (JXTJ) yield 

n ( ~0 2ig 2 io - + S a )( K 2 + u 2 ) 

C h {U},00) = n 9/ ; - ^ — ; - rr 7 | C h (u,0) 



C v (u,oo) 



2g 2 (n + iuj) - i(cu + <S q )(k 2 + uj 2 ) 



2g 2 (n + iiv) - i(u + 8 a )(n 2 + uj 2 ) 



(25a) 

C h (u,Q) 
(25b) 



In our previous paper [L4{ . we introduced coefficients 
Ck+ = (Chk + C v k)/V2, in terms of which a fidelity F 
was defined by 

e 2i *F = -2j2Ck+(0yC k+ (oo) 



= -2 J C + {uj,0)*C + {uj,oo)duj (26) 

where, in the second line, we have expressed the same 
quantity in terms of the continuous-spectrum coefficients 
we are using here. It is easy to see, using Eqs. (|25|) . that 
this expression is given by 



e 2l *F 



2# 2 (k - iuj) + i(u + 5a)(K 2 + w 2 ) 
2g 2 (n + iu) - i(w + S a )( K 2 + w 2 ) 



C h (u,0)\ 2 du 



(27) 



with 



tan 



(w + S a )(K 2 + w 2 ) - 2g 2 u 



(28) 



Clearly, in the adiabatic limit, when Ch{u, 0) is so sharply 
peaked that uj can be replaced by a constant detuning, 
A, in the phase factor in the integrand, one has F = 1 
and $ = ip(A). Additionally, in this case, as it is easy 
to see, the phase $ equals the parameter <f> appearing in 
the ideal transformation (TTJ (compare also Eq. (|28p to 
Eq. (34) of Q). 

For shorter pulses, one may write oj = A + to', ex- 
pand ([28]) around to = A, and substitute in (|27l) . Assum- 
ing a Gaussian pulse, \Ch(u>, 0)| 2 = (T/^/2~n) exp[— (uj — 
A) 2 T 2 /2] (compare Eq. (12) of [II]), after some straight- 
forward algebra, one obtains the following relatively com- 
pact form for the fidelity: 



F =1 - 



16<? 4 k 2 l-2g 2 + 2S a A + 3A 2 



T 2 (4g 4 K 2 + ((A + (5) (A 2 + k 2 ) - 2g 2 A) 



(29) 



The most interesting feature of this result is that it is pos- 
sible to make the 1/T 2 term vanish for a suitable choice 
of parameters. For instance, suppose we want to carry 



-o.u 



0.2 



0.2 

A/k 



0.1 



0.0 



FIG. 5: Behavior of the log (base 10) of the infidelity around 
the point where the second-order term in Eq. (I29[) vanishes. 
Parameters are: T = 10/ft, g 2 = (169/175)ft 2 , and 8 a as 
given by Eq. ((2}. The value of 1 — F 2 at the minimum is 
about 3 x 10" 4 . 



out a VSWAP gate. As indicated above, this means we 
want ip(A) — 7r/4, from which Eq. @ immediately fol- 
lows. Additionally, in order to make the second-order 
term in (|29[) vanish, one should set 



2g l - 3A 2 - k 2 
2A~ 



6a 



(30) 



Combining © and (|30p one gets the following equation 
for A: 

A 4 + 2A 2 ( 5 2 + k 2 ) + Ag 2 nA + k 2 (k 2 - 2g 2 ) = (31) 



It is possible to show that at least one real root to ([31 
always exists, provided 



2g 2 > (12V3 - 20)k 2 « 0.785k 2 



(32) 



that is, provided the cavity is "good enough" . 

The simplest solution to the equations (|30[) and (|3Tj) 
happens when 2g 2 = k 2 exactly. In that case, one can 
simply set A = and 5 a = k. However, as long as 
one can adjust the detunings A and S a independently, 
such precise engineering of the coupling to the cavity is 
not necessary: one only has to make sure that g is large 
enough to satisfy Eq. (pl2"j) . 

The possibility of eliminating the error to second order 
in 1/T allows for a substantial relaxation of the pulse 
length requirements. We illustrate this with the results of 
"exact" numerical calculations such as the ones reported 
in 14| . Figure 5 shows the behavior of 1 — F 2 (sometimes 
called the "infidelity") as a function of A for a pulse 
of length T = 10/k, g 2 = (169/175)k 2 , and 5 a chosen 
to satisfy ([2]) for every value of A. The figure shows 
a decrease in the infidelity of more than one order of 
magnitude at A = k/5, which is an exact root of (|3"Tj) for 
this value of g. 

Figure 6 shows the dependence of 1 — F 2 on T for 
g 2 = k and detunings chosen to satisfy Eqs. (|3"0")) and 



1.2 1.3 
log(reT) 



FIG. 6: Scaling of the infidelity with the pulse duration for 
g = ft. Both branches of solutions are shown. 



(|3"T|) . Both branches of solutions to Eq. (f3Tj) are shown. 
The best fit slopes are —3.95 ± 0.05 for branch 1 and 
-4.03 ± 0.01 for branch 2. 

In principle, for an arbitrary initial state, the "failure 
probability" of the gate (defined as 1 minus the absolute 
value squared of the overlap of the actual final state with 
the desired final state) will depend both on the deviation 
of F from unity and on the deviation of the phase $ in 
Eq. (|2"fj)) from the desired value (e.g., (j) — 7r/4 for a 
%/SWAP gate). However, it is easy to see that, since the 
failure probability has to be a positive quantity by def- 
inition, the phase error must enter it squared, to lowest 
order. Our calculations show that $ — cf> scales already 
as 1/T 2 , so that contribution to the total error is always 
at least 0(1/T 4 ). 

We can be a bit more precise. By expanding ip{ui) in 
Eq. (j2"T|) around cj = A, we get, to lowest order in 1/T 2 



1 - F 



l2 



2($ 



T 2 
y-2 " 



(33) 



where the primes denote derivatives with respect to lo. 
By the above argument, it is more important to make 
i// vanish than -0". We find that, generally, ip" does not 
vanish for the values of S a and A that cause tp' to vanish, 
except, interestingly, at the precise point where the cavity 
becomes "good enough" (according to Eq. ([32]) ). namely, 
g ~ 0.63k. This is shown in the numerical calculations 
plotted in figures 7 and 8. Except in the neighborhood 
of this point, however, it appears that the actual value 
of gj k, does not make much of a difference to either the 
fidelity or the phase error. The figures also show that, 
in general, depending on the values of T and g/n, one 
branch or the other of the solutions of (|32"|) may give a 
higher fidelity, but the differences are not very dramatic. 

Lastly, for completeness, we would like to briefly dis- 
cuss the nonadiabatic corrections to the direct SWAP 
gate. Although a SWAP between the atom and the 
photon can always be carried out via two consecutive 
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: : -x- t = io/k 

[ 1 -B- T = 5/« 

0.5 0.63 1 1.5 2 

ff/t 

FIG. 7: Logarithmic (base 10) plots of the infidelity and the 
phase error as a function of g/n, for the values of the detuning 
corresponding to one of the branches of solutions of Eq. (|31|) . 
Pulse durations of T — 10/jc, 5/k and 2/k are shown. 



VSWAP gates, it may be useful in some cases to per- 
form this operation in a single step, and, as pointed out 
in the introduction (see also [ll|) the KIN mechanism 
allows one to do that. In principle, the simplest way to 
do it is to set A = 5 a = 0, but if one also wants to elim- 
inate the second-order term in (1291) an alternative must 
be sought. As Eq. (f28|) shows, for a SWAP gate one has 
to impose (A + S a )(K 2 + A 2 ) — 2g 2 A = 0, which is satisfied 
whenever 

S ° = A (j?+~E? ~ *) (34) 

If Eq. (|34|) holds, the quadratic term in (f29|) can be made 
to vanish by the choice 



A 2 = gy/g 2 + 4k 2 - g 2 - k 2 (35) 

This equation always has a real solution provided that 
2g 2 > k 2 (the good cavity regime again). Hence, in this 
regime, high-fidelity single step SWAP gates with rela- 
tively short pulses are also possible. Numerical results 
for the infidelity are shown in Fig. 9. 




0.5 0.63 1 1.5 2 



-1 - 



-3 - 




0.5 0.63 1 1.5 2 

FIG. 8: Same as Fig. 7, for the other branch of solutions of 
Eq. ([311). 




1 1.1 1.2 1.3 1.4 

log( K T) 

FIG. 9: Scaling of the infidelity with the pulse duration for 
g = k and detunings chosen to satisfy the SWAP condition. 



VI. CONCLUSIONS 

In this paper, we have studied how to carry out cavity- 
mediated photon-photon gates by the KIN mechanism in 
a variety of systems that depart from the original pro- 
posal in a number of ways: two-sided or bidirectional cav- 
ities, nondegenerate atomic transitions, and somewhat 
longer pulses. We find that, in principle, all of these 
systems are usable, although at the expense of greater 
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experimental complexity. Alternatively, our results show 
how the KIN gates deteriorate when the original ideal 
conditions are not met, if compensating measures are not 
taken. 

A particularly interesting result is the possibility of 
choosing the detunings so as to eliminate the first correc- 
tion to the adiabatic approximation. Although one still 
needs to satisfy the condition T > 1/k, this means that 
fairly high fidelities (F 2 ~ 0.995) can be achieved for T 
as short as 5/k and improve rapidly as T increases. All 
these results increase our hope that a proof of principle 
of this scheme for photon-photon quantum logic may be 
experimentally achieved in the not too distant future. 
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Appendix: Cavity-induced phase shifts 

As mentioned at the beginning of Section IV, the 
Hamiltonian we use here (Eq. (JTTJ) ) differs from the one 
we used in [lij in a couple of particulars. One is that the 
quantization length L appears everywhere divided by 2, 
relative to the way it appeared in [l4|]; this simply re- 
flects the fact that for the traveling waves introduced in 
Section III, the natural boundary conditions are periodic, 
which lead to a mode spacing 2ttc/ L, instead of ttc/L as 
was the case for the standing waves used in [l4[. This is 
merely a formal difference, since L is, after all, arbitrary. 

A more subtle difference concerns the coupling coeffi- 
cients to the individual "modes of the universe" inside 
the cavity, which in this paper's treatment have the form 



y/2cn/L 



y/2cn/L 



(A.l) 



which differs from the form used in [14j (see, e.g., 
Eq. (11) of that paper) by a phase factor e l9k , with 
9k = tan _1 [(f2fc — Q c )/k]. This difference actually can be 
traced back to the phase choice for the incident modes in 
Section III. As suggested in Fig. 2, here we are implicitly 
assuming for the incoming field the form 



E- 



(+) 



E 



2e AL 



1/2 



whereas in [14| we had 



E: 



( + ) 



1/2 



ikL—iClkt 



(A.2) 



(A.3) 



Apart from some trivial phases, the key difference is in 
the factor £fce lfcL in Eq. (|A.3[) . By follo wing the argu- 
ments given at the end of Section III of 14|, it can be 



seen that this combination is essentially equal to e~* efc , 
with 9k defined as above. Hence, in both papers, the 
key point is that the mode operators inside the cavity 
are multiplied by a factor of e l6k , relative to the incident 
modes. We believe the phase choice made in the present 
paper is the more natural one. We should note that, in 
fact, in our previous paper our choice of an initial pulse 
spectrum (Eq. (12) of [l4|) did not properly account for 
the £,k£ lkL factor. This is relatively unimportant, since 
the choice of initial state is also to some extent arbitrary, 
but it needs to be pointed out. 

The phase choice also has consequences for the outgo- 
ing modes. In the outgoing field was written as 



E (+) _ _ 

-^out 



i v / m k \ 1/2 

2i^\e ALj 4 



-ikL—iQkt 



(A.4) 



and the point was made that the overall phase difference 
with the incoming field was a factor — e - 2lkL — —e 2 ™ k in 
our current notation. By expanding the result (|3]) for A 
in Section III above, on the other hand, we can see that 
the outgoing field in this paper can be written as 



e6 



e a k e * 



with 



(A.5) 



(A.6) 



This again shows the consistency of the two treatments 
(except for an overall minus sign which depends, basi- 
cally, on whether one treats the input mirror as having 
imaginary transmission coefficients or not), within their 
different phase conventions: the key thing to note is that 
going from the input field to the cavity field one gains 
a factor of e l8k , and going from inside the cavity to the 
output field one gains yet another one. 

These cavity-induced modifications of the spectrum 
are "built-in" in the modes of the universe formalism, so 
that, for instance, for a pulse incident on an empty cavity, 
the state vector coefficients do not change even though 
the pulse spectrum does. It can be argued, therefore, 
that the measure of fidelity that we have used here, in- 
volving only the state vector coefficients, does not really 
capture all the ways in which the interaction with the 
atom-cavity system can modify the single-photon pulses. 
This is technically true but not likely to be very impor- 
tant in practice, for several reasons. First, as long as 
the adiabatic approximation holds, the factor (|A.6[) is 
approximately independent of k, and hence in this limit 
one merely has a constant phase factor that affects all 
pulses equally upon interacting with a cavity, indepen- 
dently of the pulse polarization or the atomic state (in 
the nondegenerate case, the phase factors for vertically 
and horizontally-polarized pulses will be different if they 
have different central frequencies, but again, as long as 
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they are constant, they can be removed by a method 
analogous to the one mentioned after Eq. above). 
Second, the corrections to the adiabatic regime arising 
from the factors ()A.6|) involve the power series 



tan 



fib — fl r 



o 



(A.7) 



When this is exponentiated, the first term, being linear in 
0^, simply amounts to a displacement of the pulse in real 



space, which, being the same for all initial states (unlike 
the interaction-induced shifts discussed in Section V) can 
be easily removed or simply accounted for exactly. This 
leaves the next order terms, which are cubic corrections 
to the phase, and would end up squared as contributions 
to the error probability, so they go at least like 1/T 6 
(where, as usual, the pulse bandwidth is of the order of 
1/T). 
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